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Abstract 

The relation is studied between the entanglement production and collective radi- 
ation by an ensemble of atoms. Entanglement production is quantified by means of 
a general measure introduced earlier by the author. Primary emphasis is placed on 
the entanglement generated by pseudospin density matrices. The problem of collective 
atomic radiation can be described by the pseudospin evolution equations. These equa- 
tions define the evolutional entanglement generated by the related density matrices. 
Under conditions of superradiant emission, the entanglement production exhibits sharp 
peaks at the delay time, where the intensity of radiation is maximal. The possibility 
of regulating the occurrence of such peaks by punctuated superradiance is discussed, 
which suggests the feasibility of punctuated entanglement production. 



1 Introduction 



The concept of entanglement is at the heart of modern ideas on quantum information pro- 
cessing and quantum computing [1-5]. There exist different suggestions for manipulating 
quantum entanglement, in particular, by employing spin systems or considering atomic sys- 
tems interacting with electromagnetic fields [6]. The aim of the present paper is to analyze 
the relation between entanglement and an ensemble of radiating atoms. A specific feature 
of this analysis is the consideration of a large ensemble of many atoms under conditions of 
well developed collective effects and coherent radiation. 

First of all, it is necessary to emphasize that, generally, there are two kinds of entangle- 
ment, the entanglement of a quantum state and the entanglement produced by a quantum 
operation. And one has to distinguish between these two principally different notions. 

For pure quantum states, represented by wave functions, the concept of entanglement 
is straightforward. A pure state of a complex system is termed entangled if and only if it 
cannot be represented as a tensor product of wave functions pertaining to different parts of 
the given system. 

For mixed states, represented by statistical operators, the notion of entanglement is 
essentially more complicated. This is because functions and operators are fundamentally 
different mathematical quantities. Despite of this, the name of a state is also commonly 
applied to statistical, or density, operators. One decides whether a statistical operator (a 
mixed state) is entangled or not according to the formal structure of the given operator. An 
operator p is called separable if and only if it can be written as a convex combination of 
product states: 

p = E A ^ ®j pv ] ; X > o , a„ = i . 

V V 

When p is not separable it is named entangled. However, from the mathematical point 
of view, the properties of an operator are correctly defined not according to its formal 
appearance but respectively to its action on a given set of functions. 

In order to better understand the difference between the structure of an operator, be it a 
statistical operator or not, and its action, it is useful to distract ourselves, for a while, from 
physical applications and to clear cut the point in mathematical terms. Then, for a given 
operator, one could define two types of entanglement. One is the entanglement of an operator 
and another one is the entanglement produced by an operator. In the first case, one treats an 
operator as being entangled or not according to its formal structure of being nonseparable 
or separable. In the second case, one should characterize an operator as entangling or not 
on the basis of its action over a set of disentangled functions. Pay attention that entangled 
and entangling are quite different characteristics. 

In physical applications, one usually discusses the first type of entanglement for a partic- 
ular class of operators, that is, for statistical, or density, operators. The entanglement of a 
statistical operator is commonly termed the entanglement of state. A measure for this type 
of entanglement can be well defined only for pure bipartite systems. No general measure of 
this type is known for multipartite or mixed systems. It is worth stressing that, whatever 
entanglement of a state one is talking about, be it the entanglement of formation, entan- 
glement cost, relative entropy of entanglement, or entanglement of distillation, all of them 
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pertain to the same type characterizing the entanglement of an operator. 

The entanglement produced by an operator is a qualitatively different notion. To my 
understanding, only this type of produced entanglement, or entanglement production, can 
be correctly defined in mathematical terms. And it is just the entanglement production 
that is considered in the present paper. A general measure of entanglement production, 
realized by an arbitrary operator and for systems of any nature, be it multipartite or mixed 
systems, has recently been introduced [7,8]. This measure is studied here for characterizing 
entanglement produced in the process of collective atomic radiation. 



2 Entanglement Production 

First, we need to describe a general way of quantifying entanglement production. Consider 
a composite system consisting of parts enumerated by the index i G {i} pertaining to a label 
manifold {i}. For the standard case of a discrete manifold, one has % — 1, 2, ... In general, 
the label manifold can be continuous. Let each part be characterized by a single-partite 
Hilbert space 

Hi = Z{\ni >} , 

being a closed linear envelope of a single-partite basis {|rij >}. The composite Hilbert space 

n = ®iHi (i) 

is the tensor product, which is a closed linear envelope 

H = C{\n >} 

of a multipartite basis {|n >} consisting of the vectors 

| ra > = <g)j \rii > . (2) 

For a continuous label manifold {i} the products (1) and (2) are the continuous tensor 
products, introduced by von Neumann [9] and employed for particular cases in Refs. [10,11]. 
In each space Hi, for functions ipi and ip[, a scalar product (ipi, y^) is defined. The norm, 
generated by this product, is 

INk = \J(vi,Vi) , (3) 

which is termed the vector norm. 

The composite space (1) can be decomposed into two sets of functions. One sort of 
functions / G Ti is the class of factor functions 

f = ® m (ipi e Hi) , (4) 

which are called disentangled functions. All such functions form the disentangled set 

V = {® m \ & G Hi} . (5) 
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The remaining vectors of Ti constitute the complement Ti,\T> whose elements cannot be 
represented as tensor products of ipi G Hence, the set Ti\V can be named the entangled 
set. By construction, 

n = vun\v . 

The set V is not necessarily a subspace of Ti. Nevertheless, on the basis of the norms (3), it 
is straightforward to define the restricted vector norm 



v 



=nn^ik (6) 



over the set V C Ti. 

Let an operator A be given on 7i. For instance, this can be an operator of an observable 
quantity, or an operator representing a measurement or any other operation on H. Assume 
that the operator A is bounded, possessing the norm \ \A\ \t-c- We may introduce the restricted 
norm 



\\Af\ 



\v = sup ''| 7 M '' P , ( 7 ) 
fev 11/ 1 |d 



interpreted as the supremum of the expression corresponding to the standard definition of 
a norm over an Hilbert space, restricted to the subset V of disentangled functions. The 
operator norm \\A\\, generated by the vector norm ||A/||, is also called the Hermitian or 
spectral norm. The norm (7) can be written as 

\(f,Af')\ 

\\ A \\v = sup , (8) 

f,f'ev WfWvWrWv 



where /, /' ^ 0. For a self-adjoint operator, this can be reduced to 



|A||d = su P %^!, (A+ = A), (9) 
fev 11/ |b 



where again / ^ 0. 

Any factor function / ED can be written as an expansion 

Hi 

which can be represented as 

f = J2 c "\ n> (n = {m}), 

n 

with 

c„ = ]J a n t , \n> = <8>i \rn > 

i 

Therefore, norm (9) takes the form 

Smn Cm^n ^ 771 1 .A 1 71 > 



\A\\ v = sup . (10) 

{c„} Z^n \ C n\ 
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The latter, keeping in mind an orthonormalized basis, such that < m\n >= 5 mn , can be 
simplified to 

\t> = sup I < n|A|n > | . (11) 



This representation is the most convenient for practical calculations. 

For each operator A on Ti, one may put into correspondence a nonentangling operator 
A® having the structure of a tensor product ®iAi of single-partite operators 

A^Tr^yA. (12) 

To preserve the scale- invariant form of A®, we require the validity of the normalization 
condition 

Tr w A® = Tr n A , (13) 

in which 

Tin ® i A i = ]lTr H . A,. 

i 

The amount of entanglement produced by an operator A is characterized [7,8] by the measure 
of entanglement production, which compares the actions of A and A® on the disentangled 
set V. The nonentangling counterpart of A, taking account of condition (13), reads 

^ T Tr lf 4 (14) 
Ti n ®i A { 

The measure of entanglement production by an operator A is 



\A 



<*) = log TT^~ • (15) 



\A® 



YD 



Here the norm ||A||x> is given by Eq. (11) and the norm ||A®||x> can be calculated as 

Measure (15) will be used in what follows for particular physical applications. By construc- 
tion, this measure satisfies all properties typical of any measure [7,8]. 

Several remarks are in order. When one is interested in the entangling properties of a 
family A = {A} of operators, one may calculate measure (15) for each of them and then 
find the maximal value 

e(A) = sup e(A) 
AeA 

quantifying the entanglement production by the family A. In particular, A can represent an 
algebra of observables. 

By replacing in the above expressions the disentangled set T> by the entangled set 7i\T>, 
we can introduce the measure of disentanglement production 



S(A) = log 



\A\ 



H\D 



\\A®\\n\v ' 



describing the disentangling properties of an operator A. 

The restricted norm (11) does not necessarily coincide with the standard norm, so that, 
generally speaking, 

\\A\\ V ^ \\A\\ n , 

though in many cases they can be identical. An example, when these norms are different, 
can be composed as follows. Let ip &TI, with \ \(p\ \ = 1, can be written as 

<P = J2 c n\n > (n= {rii}) . 

n 

Define the operator 

A = \ip >< ip\ , 

for which = 1. At the same time, the restricted norm (11) gives 

H^llx, = sup |c n | 2 . 

n 

Assuming the case, when J2 n \ c n\ 2 = 1 and at least two of c n are nonzero, one has sup n \ c n \ 2 < 
1. Thence 

\\A\\ v <\\A\\ n = l. 

In general, the restricted norm | \A\ \x> can be expressed through the standard norm | \A\\n 
as follows. Introduce the projector P-p, which projects the composite space H onto the 
disentangled set V, so that Pv'H = V. Then, we may write 

\\A\\ V = \\P v AP v \\ n . 

From here, \\A\\ V < \\A\\ n . The projector P v , however, is nonlinear, hence P V AP V is a 
nonlinear operator, even if A is linear. Note that the operator norm is usually defined for 
linear operators. 

For pure bipartite systems, one measures entanglement by means of the reduced von 
Neumann entropy 

S N = -Ti Ht Pi log /3i , 

in which 

pi = Tr-H&i P , 

and p is a statistical operator of the bipartite system. Then 

S N = |Cn| 2 log|c n | 2 , 

n 

where \c n \ is a Schmidt weight. The entropy Sn is assumed to quantify the entanglement 
of a bipartite state. This should not be confused with the entanglement produced by the 
statistical operator p, which is quantified by measure (15) yielding 

e(p) = - log sup |c„| 2 . 

n 
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As is evident, these measures are, in general, different, coinciding only for the maximally 
entangled states, for which \c n \ — const. 

To stress once again that entanglement production by an operator and entanglement of 
an operator are rather different notions, let us consider the Werner operator [12], which is 
defined for bipartite systems of the composite space H — Hi ® 7^2, with the single-partite 
space dimensionality d = dim7Yj. The Werner operator is 

W a = — ^ — ; \(d - a)l + (da - l)a] , 
a 6 — a 1 J 

where a is a real number, 1 is a unity operator 

d 

1 = \niTij >< rijUil ? 

and a is a flip operator 

d 

a = \niTij >< riiUjl . 
The name of the latter operator comes from the property 

0-(</?l <g) ip 2 ) = (fi2 ® <fl ■ 

Also, one has 

Ti H W a = 1 , a = Tr w # Q( x . 

The positive partial transpose criterion [13,14] tells us that W a is separable if and only 
if a > 0. Because of this, the entanglement of the Werner operator for a > must be 
zero. However, this operator is entangling for all a ^ 1/2, which means that it transforms 
nonentangled functions f £ T> from the disentangled set V, which are the factor functions 
of the form 

d d 

f = ^Z a i\ n i > ®^2bj\nj > , 

i=l j=l 

into entangled functions W a f £ T~L \ T>. Therefore, the measure of entanglement production 
e(W a ) must be nonzero. For instance, in the case of two-dimensional single-partite spaces, 
when d — 2, we have 

1 
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From here, we find 



(2 - a)l + (2a - l)a 



W? = \l> \\W*\\v = \, \\W a \\v = l sup{2 -«,! + «} 



Then, measure (15) becomes 



:(HV) = log(l + ^|2a 



which is nonzero for all 1/2. Another example of an operator that is separable, but at 
the same time entangling, can be found in Ref. [8]. 
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3 Spin Entanglement 



The process of radiation by resonant atomic systems can be formulated in terms of the 
evolution equations for spin operators. The latter are often called pseudospin operators in 
order to emphasize that they do not correspond to actual spins but are just mathematical 
structures possessing the properties of spin operators. For resonant, effectively two-level, 
atoms, such operators correspond to spin-1/2 operators, which is implied in what follows. 

It is convenient to deal with the spin operators S? = S a (rj) being either the ladder 
operators Sf or the z-component operator S*. The former are used for describing dipole 
transitions, while the latter represents the population difference. These operators possess 
the properties 

I I 

$i $i = "F 2 ' = -'"2^ ' 

c+ o— i c z 0-0+—^ qz c± c± _ n Q z Q z — ^ 

°i J i ~ 2 * ' * * ~~ 2 ~~ ' ~~ 4 

and satisfy the commutation relations 

[Sfi Sj ] = ZdijSi , [Sf, Sj] = ^5ijSf . 
Let us introduce the spin density matrices [15] of two types. One type is 

Rn = [Rn(h ■ --injl ■ ■ ■ 3n)] , 

with the elements 

R n (i 1 ...i n ,j 1 ...j n ) = <Sl...S+S-...S-> . (16) 

And another type is 

Z n = [Z n (ii . . .i n ,ji ■ ■ -jn)] , 

with the elements 

Z n (i 1 ...i n ,j 1 ...j n ) = <Sl...S^ 1 ...Sl> . (17) 

Here the angle brackets < . . . > imply statistical averaging with a statistical operator p. Ex- 
pression (16) defines the dipole-transition, or transverse, matrix. Equation (17) determines 
the population-difference, or longitudinal, matrix. These are the matrices with respect to 
the indices i and j. 

Spin density matrices [15] can be used for describing their entanglement production under 
magnetic transitions [8]. They also appear when studying the thermal entanglement caused 
by statistical operators of magnetic systems [8,16,17]. However let us recall that in the 
present paper the spin operators will be applied for describing not magnetic systems but 
ensembles of radiating atoms. Similar pseudospin operators arise when considering atomic 
entanglement in two-component Bose-Einstein condensates [18] or in systems of atoms with 
several internal states [19,20]. 

The spatial variables rj can be either discrete or continuous, which is not essential, since 
the final results will not depend on this difference. In the case of discrete e Z 3 , one 
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employs summation over % = 1,2, ... ,N, while for continuous G M 3 , one should replace 
the sum J2i by the integral / cfav 

Single-partite spaces can be constructed on the basis of plane waves. For the discrete 
spatial variable, one has 



/AT 

while for the continuous variable, 

M^) = -j= e*"' (r< G M 3 ) , 

V" being the system volume. Since the final results do not depend on whether is discrete 
or continuous, we shall employ the simpler notation corresponding to the discrete case. The 
functions (pk(^i) are orthogonal, 

$^fc( r i)¥>fc( r j) = % > 

fc 

and form a complete basis, such that 

TV 

Hy^(r;)y?p(r;) = 8 kp . 
i=i 

Denoting the vector 

I A; >< = bfc(ri)] 

as a column with respect to i, we define the single-partite space 

n t = c{\k> t } 

as a closed linear envelope over the basis {\k >;}. Then the composite space has the form 
Ti. = <S>iHi, in agreement with Eq. (1). 

The first-order spin density matrix (16) is 

R 1 (i,j) = <S+Sr> . (18) 

For its trace, we have 

N N 

Tr^^iMM)" (l + s) , (19) 
i=i A 



where the notation 



i=i 

is introduced. 

The second-order spin density matrix (16) has the form 



2 N 



Ri(iii2jij2) = < SfSfSrSr > . (21) 
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Taking into consideration the properties of the spin operators, we see that 

[(1 - S ili2 )(l - 5 jlj2 ) - 1] R 2 (iii2jij2) = . (22) 
Studying the semi-diagonal element 

R 2 (U,jl)=R2(U,lj) , 

we shall invoke the decoupling 

< SfSfSr > = < Sf >< S+Sr > , (23) 
valid for I ^ i,j. Then we get 

R 2 (li, Ij) = (1 - 5 U )(1 - 6 jt ) Q + < Si >) < S+S- > . (24) 
This can be used for calculating the partial traces 

R{ = Tr^ 2 i? 2 , R\ = Tr^j_R 2 (25) 
of R 2 defined on Ti = TL\ <S> Ti. 2 . The elements of R®, with a = 1,2, are 

N N 

Rl(i,j) = "£R^{ilJl) , Rl(i,j) = 1^(11,1 j) ■ 
i=i i=i 

The number of atoms iV is assumed to be large, N 3> 1. We find 

N 

R a l (i,j) = -(l + s)R 1 (i,j). (26) 

From here, it is easy to get 

AT 2 

Tr n R 2 = — (1 + sf . (27) 

The nonentangling counterpart of R 2 is proportional to R\ ®R\ , with the proportionality 
constant obtained from the normalization condition 

Tin R2 = Tr w Rf . 

As a result, we come to 

flf = , (28) 

where the equality 

Ti n (R\ ® = (Ti n R 2 f 

is used and Tr n R 2 is given by Eq. (27). 

To find the entanglement produced by R 2 , we need to define the restricted norms 

||-Ri||d = sup < k\R\\k > , | |i? 2 | |d = sup < kp\R 2 \pk > . (29) 

k kp 
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For this purpose, we shall use the single-mode laser approximation for the correlator 

< S+Sr > = | e -*o*« (i ^ j) , (30) 
where k is a fixed wave vector of laser propagation and = — r^. Then we obtain 

|| J R 1 || 7 , = sup|i(l + s), . (31) 

For matrices (25) and (28), we find 

ll-^illf = ~2 s )ll-^illf > ll-^f I |x> = ll-Ri|lx> • (32) 
The measure of entanglement produced by R 2 reads as 

e(R 2 ) = log . (33) 



To work out the norm H-R2II1), we shall involve the decoupling 

< Sj~ S m S n > = < S^S n >< Sj~S m > + < S^S m >< Sj~S n > , (34) 
valid under condition that all indices are different. Then we get 

|| J R 2 || 2 > = sup|| (1 + s) 2 , ^w 2 | . (35) 

Comparing Eqs. (31) and (32), we see that 

I I-R2I \v — 2| I© • (36) 

Therefore norm (33) reduces to 

e(R 2 ) = log 2. (37) 

Thus, the entanglement produced by the transition spin density matrix R 2 is always constant. 
Now, let us turn to the longitudinal spin density matrices (17). Such a first-order matrix 

is 

Z 1 (i,j) = <S]S?> . (38) 

As is evident, 

N N 
i=i 4 

The second-order longitudinal spin density matrix writes as 

z 2 (t 1 t 2 ,j l32 ) = <s* 2 s; i s* i s* 2 > . (40) 

The matrix Z 2 is defined on H = Hi <g> H 2 . The partial traces give the matrices 

Z\ = Tr m Z 2 , Z\ = Tr m Z 2 , (41) 
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whose elements are 

N N 

Zl(i,j) = ]TZ 2 (^7) , Z^iJ) = £Z 2 (K,y) ■ (42) 
1=1 1=1 

We find 

Z 2 (li, Ij) = 1 < 5/5? > , j) = ^ < 5/5? > , (43) 

for any a = 1,2. From here, 

TV 2 

Tr H Z 2 = — . (44) 
lb 



The nonentangling counterpart of Z 2 is 



For the restricted norms, we have 

\\Zi\\v = \ (1 + iVs 2 ) , ||Z®||x> = \\Z x \f v . (46) 
The measure of entanglement produced by Z 2 is 

e(Z 2 ) = log . (47) 

Substituting to Eq. (47) the norm 

||Z 2 || E = SU p{i_,^}, (48) 

we see that the entanglement production (47) strongly depends on the value of s, which is 
the average population difference. When the latter is zero, 

e(Z 2 ) = \og3 (s = 0). (49) 

But for any nonzero s, we have 

e(Z 2 ) = (s ± 0) , (50) 

in view of large N >> 1. The temporal behaviour of s = s(t) depends on particular physical 
systems. 



4 Collective Radiation 

Our aim is to find the temporal evolution of the average population difference s(t) under 
collective radiation of atoms. In this section, the derivation of the pseudospin equations 
governing the behaviour of s(t) is given. This is done for resonant, effectively two-level, atoms 
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with dipole radiation transitions. One usually assumes that electrodipole and magnetodipole 
transitions are described by equations possessing identical structure. However there exists a 
difference between the equations for these two types of transitions. In order to analyze this 
difference, we shall consider here the case, when electrodipole transition is forbidden, and 
the radiation occurs at magnetodipole transition. 

The Hamiltonian of an atomic system, emitting electromagnetic radiation, is 

H = H a + H f + H af . (51) 

The first term 

N /I \ 

(52) 



corresponds to resonant atoms, with transition frequency uj . The field Hamiltonian writes 

as 

Hf = ± J (E 2 + H 2 ) dr , (53) 

where E is the electric field and H = V x A is the magnetic field. The vector potential A 
satisfies the Coulomb calibration 

V • A = . (54) 

The last term H a f describes atom- field interactions [21]. 

When electric dipole transitions are forbidden, we need to consider other radiation tran- 
sitions. We study here the case of magnetic dipole transitions. The matrix element of the 
magnetic dipole 

with the light velocity, c, is expressed through the density of current 

jmn(r) = - (CV^n ~ ^Vft) » 

in which e is electric charge, m e its mass, and ip n = ip n (r) are wave functions of states 
enumerated by the index n— 1,2. Electric dipole transitions are forbidden under conserved 
parity, when 

C(-r)lM-r) = COOlMr) (56) 
for all m and n. Thence, the density of current is antisymmetric, 

jmn(-r) = -jmn(r) • (57) 

Denote the transition magnetic dipole 

A* = A*2i , AT = A*i2 • (58) 

Diagonal elements of dipole (55) are also nonzero, jl nn ^ 0. This is contrary to the case of 
electric dipole transitions, when parity is not conserved, and only the nondiagonal dipole 
elements survive, with diagonal ones being identically zero. The existence of nonzero diagonal 



13 



elements ft„ is in the basis of the differences between the following equations for electric 
and magnetic dipole-transition radiation. 

Under magnetic dipole transitions, the atom-field interaction Hamiltonian writes as 

N 

H af = -Y, M *- B * + Kf, (59) 
1=1 

where the operator of magnetic moment is 

M i = /rs,+ + /rsr, (60) 

and the last term is 

= - E (ft + A/I S?) • B 4 , (61) 

i=i 

with the notation 

ft = 2 ^ n + ^ ' = ^ 22 ~~ ^ n • ( 62 ) 

The total magnetic field is the sum 

B, = H + H, (63) 

of an external magnetic field Ho and the radiation field H« = H(rj,t). 

Using the commutation relations for the operators E, A and H, which can be found in 
book [22], we have the Heisenberg equations 

1 <9E An . 1 OA 

in which the density of current is 

N 

j = -c E (Mi + ft + A/2 5?) x V5(r - r,) . (65) 
i=i 

The Heisenberg equations for the pseudospin operators are 

= -iSf (w„ - A/2 ■ Bi) - 2<5?/I • B* , = i (ft?+ - jTSr) ■ B, . (66) 

One usually considers the system of equations (64) and (66) for the field and pseudospin 
operators. However, it is possible to eliminate the field variables and to reduce the problem 
to the equations for only pseudospin variables [21]. To this end, from Eqs. (64) we get the 
equation 



whose solution reads as 



A(r, t) = A mc (r, t) + - [ j f r', t - - — — ] , (6S 



c 
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where A vac is the vacuum vector potential. Retardation in the dependence of spin operators 
can be taken in the Born approximation 

S~ (t- S-(t)Q(ct - r) e tk ° r , S* (t - ^ S*(t)e(ct - r) , (69) 

where 0(-) is the unit-step function. 

Substituting the density of current (65) into potential (68), we have 

A, = At + A" + A[ + A vac , (70) 

(71) 



where 



. r~. v c 



A 'M = ~ E "T r « X ^ + A ^ 5 I (* 



with the notation 
Here the equality 



° S7 (t- -)= ik S~ (t 



dr 3 V c/ J V c, 

is taken into account, resulting from the Born approximation (69). 
For the magnetic field 

H,(t)=H(r,,t) = V,x A t (t), (72) 

we obtain the expression 

H, = H+ + H- + + H mc , (73) 

in which 



H+W = E 



^ jl- (/x-nijQnij _ ^ + ^ ^ /I - 3(/Z • n, ; !ii; ; 



r './ 'f> 



( t - ^ 



H , ((H _ £ M--3(M-.n, jK | (?4) 



where = r^/r^- and 

r. 



M; = /r + A/2 5| [t- If) . (75) 
In what follows, we shall replace summation by integration according to the rule 



N 

E p / dr ' 



in which p = N/V is atomic density. Integrating over the spherical angle Jl(n), related to 
the unit vector n = r/r, we shall employ the properties 

8tt 



/ n)n] dfi(n) = y /Z , 
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J [/7 — 3(/7 ■ n)n] d{l(n) = , ■ n)n d{l(n) = ^ ft 



(76) 



In sums (71) and (74), there are the terms, where j = i and = 0, which yields 
divergence. These are the so-called self-action terms, whose divergence occurs because of 
the point-like representation of atoms. To avoid this divergence, the self-action must be 
treated more accurately, which can be done as follows. 

Consider a single atom, located at r.,- = 0. At the point r, it creates the vector potential 



A s — A* + A s + A s , 



in which, according to Eqs. (71), 



A+M) = -l±^rx^+(t- - 



A' 8 (r,t) = - — r x 



/7 + A/I S z ( t 



The corresponding magnetic field 

h s eVxa s = h: + h; + h;, 

contains, in agreement with Eqs. (74), the terms 



H s + M) 



kl tliE^h _ (1 + ikor) ^-3(^-n)n 



H' s (r,t) = - 



M' - 3(M' • n)n 



where 



M' = /7 + A/7 S z (t - 



To consider the self-action, we need to go to the limit r — > 0. Before doing this, let us 
average H s over spherical angles, which gives 



flS + (t- - ) + ITS- ( t 



H s ^- Jll s (r,t) dn(r) = -^ 
Keeping in mind r — > 0, we may present the exponential e* fe ° r as an expansion 

e ik r _ x + ikor ( r _> 0) . 

Imitating the atom nonlocality, we may average the singular term 1/r over the distance r 
between the electron wavelength A e = 2nTi/m e c and the radiation wavelength Ao = 2tic/ujo, 
so that 

1 



An — A fi J X 



x o dr ho I m P c 
— = — In 
r 2ir 



where we take into account that A e <C A . Then, for the averaged in this way self-action 
magnetic field, we find 



H s (t) 



fTS-(t) +jlS + (t) 
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Introduce the notation for the natural width 



and the Lamb frequency shift 



7o = \ Wh% (77) 



? h (^) ■ (78) 



Equations (66), for a single atom with no external fields, are transformed into the system of 
equations 

^1 = _ Sl _ i lo)S - - ? ( 7o + iS L ) S + + ^-l ( 7o + i8 L ) (\ - S> 
at |/z| z |/i| z \2 

dS 



= "270 (V + i) . (79) 



dt 

These equations show that the influence of the radiation self-action reduces to the appearance 
of the longitudinal attenuation 71 = 270 and transverse attenuation 72 = 70 When the atom 
is immersed into a medium, one usually treats 71 and 72 as independent parameters, which is 
assumed in what follows. The Lamb shift can be included in the definition of the transition 
frequency ujq. 

Thus, the influence of self-action will be taken into consideration by incorporating the 
related attenuation terms with 7 X and 72 into the evolution equations. And in the magnetic 
field H = V x A, the self-action will be excluded. Then the field H = H(r,t), in view of 
Eq. (72), can be rearranged to 

H = H„ on + Hdip + H mc , (80) 

where the nondipole part 

H rao?1 H non + H non (81) 
is due to the field created by radiating atoms, without the dipole part 

H, IP = H+ p + + H' . (82) 

In the corresponding expressions for the nondipole field, 

H+„(r, t) = £ kl /2 ~ (/2 - n " Kj Sf (t - (83) 
and for the dipole part 

the self-action terms are omitted. The last term in Eq. (82), H'(r, t), for r = r^, is given by 
H'i(t) = H'(rj,t) from Eq. (74), again with no self-action terms. 
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Introduce the radiating symmetric part of the total radiation field, 

H rad = ^J (H no „ + H dip ) dQ{r) = ±-j H„ on (r, t) dtt(r) . (85) 
This can be decomposed as 

H rad = H+ d + H~ ad , (86) 

with 

H r + arf (r 4 , t) = \ klfl £ -L S+ (t - T -f) . (87) 

Hence, the total radiation field can be represented as 

H = Hrad + AH + H mc , (88) 

with the notation 

AH = H no „ + H dip — H r ad ■ (89) 
The anisotropic part of the radiation field contributes to the expression 

£(r, t) = -2ifi ■ (H mc + H„ on + B_ dip - H rad ) , (90) 

for which 

Jt(r,t) cM(r) = 0. 

Equation (90) describes local field fluctuations. We also define 

Ur,t) = -Afl-H(r,t) , (91) 

where H is given by Eq. (88), and the frequency shift 

A = -A/7 • H . (92) 

Then we consider the statistical averaging of the pseudospin equations (66) in the frame 
of the scale-separation approach [21]. In this course, we define the statistical averages for 
the transition function 

u(r,t) = 2 < S~(r,t) > , (93) 

coherence intensity 

w(r,t) = A<S + (r,t)S-(r + 0,t>, (94) 

and population difference 

s(r,t) = 2 < S z (r,t) > , (95) 

in which the angle brackets < . . . > imply the averaging over the pseudospin degrees of 
freedom, not involving the fluctuating fields (90) and (91) treated as random variables. 
Function (94) is to be understood as the limit 

w(r,t) = 4 lim < S + (r, t)S~(r', t) > . 

|r— r'|— »+0 
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In the mean-field approximation, the latter simplifies to w — \u\ 2 . Definitions (93), (94), 
and (95) are in agreement with Eqs. (20) and (30). 
Also, introduce the effective field acting on atoms as 

f = fo + f rad + Z, (96) 

in which 

f = -2if2 -H (97) 
is due to an external magnetic field; the term 

f rad = -2i </7-H m(i > (98) 

is caused by the isotropic radiation field (85); and £ is the local fluctuating field (90). 
Allowing for Eqs. (86) and (87), field (98) takes the form 



frad(r,t) = -ry p J 



-»2 

G(r - r', t)u(r', t ) + ^-G*(r- r', t)u*(r', t) 



dr' , (99) 



with the transfer function 

exp(ifc r) 
= k^r ~~ ' 

Finally, from Eqs. (66) we derive the evolution equations 

9u ., . . . . 
— = -i{u + A + £ - vy 2 )u + fs , 

^ = - 2l2 w + (u*f + f*u) s , 

§ = -^7 + r«)-7i(s-C), (100) 

where ( is a stationary population difference of an atom. These equations differ from anal- 
ogous equations for electric dipole transitions [21,23] by the presence in Eqs. (100) of the 
frequency shifts A and £ - Both these shifts are due to A/7, given in Eqs. (62), which 
is nonzero because of the nonvanishing diagonal elements /7 nn . The constant part of the 
summary shift A + £o can be neglected as compared to u . And the fluctuating part of this 
shift results in the appearance of an additional inhomogeneous broadening 73, so that the 
total line width becomes the sum 72 + 72- 



5 Entanglement by Superradiance 

According to Sec. 3, the entanglement-production measure e(Z 2 ), caused by the spin density 
matrix Z 2 , defined in Eq. (17), strongly depends on the value s = s(t) of the population 
difference. Changing s(t) would govern the level of produced entanglement. To make such 
an evolutional entanglement efficient, one needs that s(t) would cross the value s = 0, when 
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the measure s{Z^) jumps between expressions (49) and (50). Thus, we have to analyze the 
temporal behaviour of the population difference s{t), which is described by Eqs. (fOO). 

It seems that an adequate method of quickly changing the population difference between 
s 7^ and s = could be by means of the process of superradiance. Then, at the initial time 
t — 0, one prepares the atomic system in the inverted state, with s ~ 1. The superradiance 
burst is peaked at the delay time to, when s ~ 0, after which s(t) goes rapidly to s = — 1. 

Temporal dynamics of superradiance, described by Eqs. (100), has been investigated for 
atomic systems [21,23,24] as well as for spin superradiance [21,25,26]. Therefore, we shall 
not repeat here the solution of Eqs. (100) but will straightforwardly go to results. 

Let the atomic system at the initial time be inverted, with s = s(0) > 0. The effective 
coupling parameter of atomic interactions through radiation field is 

7o f sm(k r - kz) 
9 = 9— / 7 dr , 101 

72 J k r 

where k = uj/c is the wave vector of the seed field selecting the longitudinal propagating 
mode. The frequency u of the seed field is in resonance with the atomic transition frequency 
uo- Because of u uo, one has k ~ k . Substantial coherence develops in the system if the 
coupling parameter (101) is large, such that gs ^> 1. 

Consider a purely self-organized process, when no coherence is imposed upon the system 
at t — 0, so that w = w(0) = 0. Then the relaxation process starts with the incoherent 
quantum stage of spontaneous radiation, which lasts till the crossover time 

t c = — , (102) 

where T2 = 1/72- For gso ^> 1, the crossover time is small, t c T2. After time (102), 
coherent effects become important and a superradiance pulse arises, with the intensity of 
radiation peaking at the delay time 



t = t c 1 + In 



l3tc 

Here 73 is the dynamic inhomogeneous broadening defined as 



(103) 



73 = Re lim - f dt /"* « f*(f)f(0 » ex P {-fao + I2 - l2gs){t - t')} dt' , 

where <C . . . ^> implies the stochastic averaging over the random fluctuating fields (90). At 
the delay time (103), the value 7 3 t c , approximately, is 

1 

73* c ~ < 1 > 

2^50 

because of which to > t c . The superradiant burst is rather narrow, with the pulse time 

r p = — . (104) 

gs 
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At the transient coherence stage, when t c < t < Ti, where 7\ = I/71, the population 
difference behaves as 



The first term here is small, provided that g ^> f . After the time to, function (f05) rapidly 
diminishes to the value about s = —s . The population difference passes through s = at 
the time, which is very close to to- Thus, the superradiant regime makes it possible to realize 
a very rapid passage of the population difference through zero. 

Then one can invert the system again by a 7r-pulse. In the process of this inversion, 
s also crosses zero. Upon regaining inversion, a new superradiant regime develops, and s 
crosses zero after the delay time t , counted from the moment of the novel inversion. This 
procedure can be repeated many times. Each time, when s = 0, entanglement production 
reaches its maximum, with the measure £(-£2) jumping from zero to log 3, in accordance with 
Eqs. (49) and (50). These jumps are so sharp because of the large number of atoms A 1 
and, respectively, owing to well developed coherent effects. A series of superradiant bursts, 
generated by the repeated preparation of an inverted atomic system, can be called punctuated 
superradiance. Such a regime can also be realized for the case of spin systems, resulting in 
punctuated spin superradiance [26,27]. This regime can be employed for regulating the 
related entanglement production and for generating a punctuated series of sharp signals of 
produced entanglement. The described regime allowing for the generation of a regulated 
series of sharp entanglement pulses can be called the punctuated entanglement production. 

The nontrivial behaviour of evolutional entanglement, considered in this paper, is due 
to the entanglement production e(Z 2 ) caused by the second-order spin density matrix Z 2 . 
Generally speaking, we could study as well the entanglement produced by higher-order spin 
density matrices Z n . However, emphasis has been placed on Z 2 because of the special role 
of second-order density matrices in statistical mechanics [28] and for quantum information 
processing [29,30]. 
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